Voltage controlled spin precession in InAs quantum wells 
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In this work we investigate spin diffusion in InAs quantum wells with the Rashba spin-orbit 
coupling modulated by a gate voltage. The gate voltage dependence of the spin diffusion under 
different temperatures is studied with all the scattering explicitly included. Our result partially 
supports the claim of the realization of the Datta-Das spin-injected field effect transistor by Koo et 
al. [Science 325, 1515 (2009)]. We also show that the scattering plays an important role in spin 
diffusion in such a system. 
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In the past decades, a great deal of effort has been 
made for the realization of the spintronic devices The 
spin-injected field effect transistor (SIFET), proposed by 
Datta and Das in 1990, 5 is one of the most intrigu- 
ing devices^— but posts some challenges to experiments 
(e.g., the spin-polarized injection and detection). Very 
recently Koo et al£ reported that, in InAs quantum 
wells with nonlocal spin valve configuration, the non- 
local voltage was observed to oscillate with the varia- 
tion of gate voltage at the low temperature when the 
two ferromagnetic electrodes (spin injector and detec- 
tor) are magnetized along the spin diffusion direction. 
They claimed that they have realized the SIFET be- 
cause the oscillation can be fitted by a theoretical equa- 
tion describing the SIFET. Nevertheless, as pointed out 
by Bandyopadhyay^ the theoretical equation adopted 
by Koo et al. only applies to the one-dimensional sys- 
tem instead of the two-dimensional one. Therefore the 
agreement between this equation and the experimental 
data 9 - makes little meaning and doubt is cast on the con- 
clusion presented by Koo et al.. Later, Zainuddin et 
aZ^i extended the one-dimensional theory to the two- 
dimensional case with an equation similar to the one 
obtained from the one-dimensional theory. However, as 
further reported by Agnihotri and Bandyopadhyayji 2 - the 
experimental data actually do not match the equation for 
the two-dimension SIFET. Therefore, whether the de- 
vice proposed by Koo et al. realizes the SIFET is still 
under debate. It is noted that all the theoretical works 
mentioned above were performed without any scattering. 
However, the scattering exists in reality and can be very 
important for spin diffusion^ ' 13 ' 14 



In fact, a thorough understanding of spin diffusion in 
the two-dimensional SIFET with the scattering explic- 
itly included can be obtained based on the kinetic spin 
Bloch equation (KSBE) approach? 4 - which has been suc- 
cessfully applied to study the spin diffusion/transport 
in various two-dimensional systems (e.g., GaAs quantum 
well o 13 ' 15 " — and Si/SiGe quantum wellsi 4 -). In the frame- 
work of this approach, spins of electrons with wave- vector 



k precess in spatial domain with frequency 

Wk = m*(fi k + gHB&)/k x (1) 

during the spin diffusion j 13 ' 15 Here, the spin diffusion di- 
rection is set to be the x-axis, m* is the effective electron 
mass, f2k is the D'yakonov-Perel' (DP)H spin-orbit cou- 
pling term and B is the external magnetic field. In InAs 
quantum wells, the Rashba spin-orbit coupling^ dom- 
inates and thus f2k = 2cv(— k y , k x , 0) with a being the 
Rashba coefficient modulated by the gate voltage. More- 
over, the small external magnetic field used to magnetize 
the electrodes can be neglected when compared to the 
Rashba spin-orbit coupling. 9 Therefore, the spatial spin 
precession frequency 

w k = 2am* (-tan k , 1,0) (2) 

depends on the polar angle #k of the momentum. This k 
dependence of the precession frequency leads to the inho- 
mogeneous broadening i 15 ' 20 The inhomogeneous broad- 
ening itself causes reversible spin relaxation during spin 
diffusion.— One notices that the one-dimension model 
adopted by Koo et al£ actually excludes the inhomo- 
geneous broadening by neglecting the transverse compo- 
nent of the momentum (i.e., k y ) and therefore is inap- 
propriate. The scattering also plays an important role 
in spin diffusion which makes the relaxation irreversible 
and affects the spin diffusion length or even the precession 
frequencyi 13 ' 14 Furthermore, the temperature (T) depen- 
dence of the spin diffusion length should be mainly from 
the temperature dependence of the scattering in this case 
as the inhomogeneous broadening is insensitive on T (cjk 
does not depend on the magnitude of k). In this paper, 
we numerically solve the KSBEs under the DP mech- 
anism and obtain the gate-voltage dependence of spin 
polarization at the detection point. Our result is in good 
agreement with the experiment of Koo et al. at low and 
sufficiently high temperatures and hence is in favor of 
their claim of the realization of the SIFET. Moreover, 
the role played by the scattering in spin diffusion is also 
investigated and revealed to be important. 

We start our investigation from InAs quantum wells 
as presented in Ref.^. The depth Vq and width a of 
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the square well are set to be 430 meV and 2 nm, re- 
spectively. The initial spatially uniform electron den- 
sity N e is 2.5xl0 12 cm -2 and the effective electron mass 
m* = 0.05too where mo is the free electron mass. The 
x-axis polarized spins (the polarization Pq is set to be 
0.02) are injected at the left boundary x — and diffuse 
along the x-axis. Due to the narrow well width, moder- 
ate electron density and small polarization, only the low- 
est subband is relevant in our investigation. The Rashba 
spin-orbit coupling coefficient a is taken from Koo et al..— 
The impurity density Ni is estimated to be 0.0397V e ac- 
cording to the mobility /i = 50000 cm 2 V -1 s _1 reported 
by Koo et al. It is noted that these parameters (e.g., 
Ni and N e , etc.) are obtained from the low temperature 
case of the experiment of Koo et al.. Although these pa- 
rameters may vary with the temperature^ we still apply 
them to the high temperature investigations due to the 
lack of necessary experimental information. The other 
parameters can be found in Ref.l22l The KSBEs read 
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Here, p k (a;,t) are the single-particle density matrices 
of electrons with the in-plane wave-vector k at posi- 
tion x and time t. *f?(x,t) is the electric potential sat- 
isfying the Poisson equation V 2 *(x,t) = e[n(x,t) — 
No]/(ano£o) with n(x,t) — X)kT r [(°k(£, *)] standing for 
the electron density at position x and time t, Nq the 
background positive charge density, and Kq the rela- 
tive static dielectric constant. — i[f2 k ■ er/2, pu.(x, t)] 
is the coherent term describing the spin precession. 
9 Qt'^ \ scg + is the scattering term with the electron- 
impurity, electron-acoustic/longitudinal optical phonon, 
and electron-electron scatterings included. The details 
of the scattering term can be found in Refs.0, [H and 
l24l . It is noted that no fitting parameter is needed in our 
calculation. 

To solve the KSBEs, the initial conditions are set as 

p k (0, 0) = (F£ t + F k °J/2 + (F k ° t - (4) 
p k (x>Q,0) = (F£ f + F£ l )/2, (5) 



and the boundary conditions are given ai 
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Pk(0,t)k> = (F£ t + F k °J/2 + (f* - ^K/2,(6) 
p k (L,t)\ k „ <Q = (F£ t + F£ i )/2, (7) 
*(0,t) = *(£,<) = 0. (8) 

Here, x = L stands for the right boundary with L much 
longer than the spin diffusion length. F k j L (F^) stand 
for the Fermi distributions of electrons with spin paral- 
lel (antiparallel) to the x-axis determined by the tem- 
perature and the initial polarization at the two bound- 
aries. The numerical scheme for solving the KSBEs can 
be found in detail in Ref.dU With the single-particle 



density matrices obtained by solving the KSBEs, the spin 
polarization at the point x at the steady state can be ob- 
tained as 

P(x, +00) = ^ Tr[p k (a;, +oo)a x ]/n(x, +00) 

k 

= ^F k (x,+co). (9) 



Since the nonlocal voltage measured in the experiment 
is proportional to the spin polarization at the detection 



point, 



we fit the experimental data measured at Xq 



with P(xq, +00). 

In Fig. [TJ we plot the gate voltage Vq dependence of 
the spin polarization at the detection point xq — 1.25 (im 
by the solid curves and that of the experimentally mea- 
sured nonlocal voltage by the dashed curves under differ- 
ent temperatures. It is noted that the Rashba spin-orbit 
coupling coefficient used in our calculation with the gate 
voltage from —5 V to —3 V is obtained by linearly ex- 
tending the a-Va curve presented by Koo et al.£ From 
Fig- HI one finds that our result is in good agreement with 
the experiment at 7 K and 40 K. Moreover, more than 
one period of oscillation are found in our result, suggest- 
ing that Koo et al. may also observe more periods if 
they enlarge the scope of measurement. When the tem- 
perature is higher, more impurities will be ionized and 
both the electron and impurity densities will increase. 21 
This may explain the discrepancy between our theoreti- 
cal result and the experimental data at T — 77 K. When 
T = 300 K and hence the scattering is strong, both our 
calculation and the experiment by Koo et al. show that 
the oscillation of the spin polarization/nonlocal voltage 
disappears. In fact, due to the suppression on spin dif- 
fusion (mainly caused by the strengthened scattering, as 
revealed in the following), the spin diffusion length be- 
comes much shorter than the spacing between the injec- 
tor and detector, and therefore no spin polarized signal 
can be observed at the detection point. 

To further investigate the influence of scattering on 
the spin diffusion, we first consider the much simplified 
case without the scattering and electric field, the single- 
particle density matrix for any k in the steady state can 
be obtained easily from the KSBEs asi£ 



— *p k (0,0)e— ^ 



PkOM), 



> 
< 



,(10) 



where u; k is given in Eq. ([2]). Then at the detection point, 



P k (x ,+oo) = 



.0, 



+ (1-* 2 ) cos(^)], fcx>0 



k x <0 



with s — ky/k — sin# k , 6 Xa — 1m* ax^ and J3 k = 
(F^ — F^)/N e . This solution with k x > has the same 
form as the result from Zainuddin et al. [Eq. (5a) in 
Ref.|Tl|. Our result clearly indicates that the contribu- 
tion to the total spin-polarized signal mainly comes from 
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FIG. 1: (Color online) Gate voltage dependence of the spin 
polarization obtained from the KSBEs (solid curves with the 
scale on the right hand side of the frame) and that of the non- 
local voltage measured in the experiment 9 (dashed curves) at 
the detection point xo = 1.25 pm under different tempera- 
tures. The plots are shifted for clarity as Koo et aZ..— 



the fc^-positive states around the Fermi circle. Instead of 
summing Pk(^o I +°°) over the fe^-positive Fermi circle 
line as done by Zainuddin et aL,— we take into account 
all the fer-positivc states and obtain 



P(xq, +oo) oc 



d0k(l 



2 m ax 2 

2 sin — cos ok)- 12 

cos &k 



It is noted that the integration over 9^ in Eq. (IT21 stands 
for the interference among different k states. However, 
one finds that this equation can not fit the experimental 
data very well as the situation faced by Agnihotri and 
Bandyopadhyay 12 until the scattering is included as pre- 
sented previously. 

We then solve the KSBEs numerically by varying the 
impurity density artificially. Without losing generality, 
we take the temperature to be 7 K and the gate volt- 
age to be 0. Under these conditions, the x dependence of 
the spin polarization with different impurity densities are 
plotted in Fig. [2] It is noted that the curve with = 
in this figure is obtained directly from Eq. (fl~2"|) . The 
impurity density dependence of the spin diffusion length 
is plotted in the inset of Fig. [5] From this inset, one 
finds that the spin diffusion length decreases sensitively 
with the increase in the impurity density. It is noted 
that even for the case of N t = 0.5N e , the system is still 



in the weak scattering limit as ojlt p — 1.19 > 1, where 
ujl = 2akp is the spin precession frequency due to the 
spin-orbit coupling and t p is the momentum relaxation 
time. The decrease in the spin-diffusion length with the 
increase in the impurity density can be understood al- 
ternatively by means of the quasi-independent electron 
model)21~— where the spin diffusion length is character- 
ized by y/DsTs with t s standing for the spin relaxation 
time and D s representing the spin diffusion constant. 
D s decreases with the increasing scattering strength^ r s 



has the same tendency as D s 

2 



as long as electrons are in 




FIG. 2: (Color online) x dependence of P with different im- 
purity densities. The impurity density dependence of the 
spin diffusion length is also plotted in the inset. The case 
of Ni — 0.0397Ve corresponds to the experimental situation. 
The temperature is 7 K and the gate voltage is 0. 



the weak scattering limits Therefore, the spin-diffusion 
length decreases with the increase in scattering strength. 
This explains the disappearance of the oscillation of the 
spin-polarized signal at T = 300 K in Fig. [T] since the 
electron-phonon scattering is strengthened there. 

In summary, we have investigated the spin diffusion in 
n-type InAs quantum wells with the scattering explicitly 
included under the DP mechanism. The consistency be- 
tween our theoretical result and the experimental data 
partially supports the claim by Koo et a/.— that a SIFET 
has been demonstrated. The essential role played by the 
scattering is also revealed. It is shown that the spin dif- 
fusion length decreases with the increase in the impurity 
density in the weak scattering limit. 
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